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Abstract: We study the effect of nonlocal optical response on the 
optical properties of metallic nanowires, by numerically implementing the 
hydrodynamical Drude model for arbitrary nanowire geometries. We first 
demonstrate the accuracy of our frequency-domain finite-element imple- 
mentation by benchmarking it in a wide frequency range against analytical 
results for the extinction cross section of a cylindrical plasmonic nanowire. 
Our main results concern more complex geometries, namely cylindrical and 
bow-tie nanowire dimers that can strongly enhance optical fields. For both 
types of dimers we find that nonlocal response can strongly affect both the 
field enhancement in between the dimers and their respective extinction 
cross sections. In particular, we give examples of maximal field enhance- 
ments near hybridized plasmonic dimer resonances that are still large but 
nearly two times smaller than in the usual local-response description. At the 
same time, for a fixed frequency the field enhancement and cross section 
can also be significantly more enhanced in the nonlocal-response model. 
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1. Introduction 



In plasmonics |[T1|2]|3], sub wave length metal structures are used to confine and enhance HHJl, 
guide ||6]|7l, and scatter |I8]|9][T0][TT] light. The goal to measure and control light at ever smaller 



length scales drives the research field towards true nanoplasmonics. It is then a natural question 
to ask: down to which sizes is the optical description of the metal solely in terms of its bulk 
dielectric function still accurate? 

One phenomenon beyond this usual description that becomes important is nonlocal optical 
response IT] [l2l El [H [15] : the fact that not only light but also moving electrons in the metal 
transport energy. Here we will focus on effects due to nonlocal response, using the linearized 
hydrodynamic Drude model HI [16] [iTl [l4l . Compared to the usual local-response Drude the- 
ory for free electrons, this hydrodynamic theory has the Fermi velocity of the electrons as an 
additional parameter 

Nowadays, simulating local-response plasmonic properties has almost become a standard 
task, even for complex geometries, thanks to the availability of advanced numerical methods 
and dedicated software. Because of the continuing success in fabricating ever smaller plas- 
monic nanostructures, the nonlocal response will become increasingly important. This has 
stimulated us to develop a similar reliable and easy-to-use numerical tool also for nonlocal 
response, and apply it to geometries where we expect effects of nonlocal response to be sig- 
nificant, as presented here. The nonlocal calculations are numerically more challenging, since 
the Fermi wavelength of the order of ' 1 nm enters as a new length scale of longitudinal 
waves HI [16] [17] [l4l . The computational grid with typical separations Ax should resolve not 
just the sub-wavelength features, but rather sub-Fermi-wavelength features of the geometry of 
typical size L and field distributions for optical wavelengths A. The numerical grid size Ax must 
be smaller than all physical length scales in our study, and the latter satisfy fcp ' <L<^X. 

The core of this article is a numerical study of nonlocal-response effects when light scatters 
off nanoplasmonic dimer structures, which are archetypical structures to study both field en- 
hancement li4l [r8ll . scattering HJIH, and hybridization of plasmonic resonances |[T9ll20l 12111221 . 
Important is also that dimers can display resonances in the visible li23]| even when their two con- 
stituents, taken separately, would not. As is known from local-response hybridization theory, 
hybridization energies grow as dimer distances are reduced. Here we study how nonlocal re- 
sponse affects his behavior. Our study of nanowire dimers complements recent work on dimers 
of nanospheres llT2l |24] [TSl . Here we present results for dimer separations only down to 1 nm, 
because for smaller separations, quantum effects not taken account into our model are predicted 
to strongly reduce hybridization energies 1251 . We focus solely on extinction cross sections and 
field enhancements that can be probed with light, leaving for later study the dark modes of the 
nanowire dimers that could be seen in electron energy loss spectroscopy 12311261 or cathodolu- 
minescence experiments l27l . 

The structure of this article is as follows. In Sec.|2]we introduce the theoretical formalism, 
and the numerical implementation in Sec. [3] In SecHwe benchmark our implementation of the 
hydrodynamical Drude model against the analytically solvable problem of the scattering off a 
single cylindrical nanowire, and also study the size dependence of nonlocal effects. In Sec.|5]we 
compare nonlocal response against local response for dimers consisting of two such cylindrical 
nanowires where field enhancement occurs in the open cavity between the cylinders. Then in 
Sec.|6lwe do the same for bow-tie dimer nanowires, where field enhancement occurs near the 
almost touching sharp tips of the triangles. Our conclusions are given in Sec.|2] and details on 
our numerical calculations in the Appendix. 

2. Theoretical formalism 

The usual local-response dielectric function of realistic metals is the sum of a Drude free- 
electron response plus interband effects IS [TTl . In the hydrodynamic description, only the 
free-electron response which is modified while the other (interband) effects are unaltered. In 
the Maxwell wave equation for the electric field E, the free-electron response is described by 



a current density J, while the rest of the optical response is modeled with a local, usually spa- 
tially piecewise constant, dielectric function eothei (r, (o). We are interested in the linear optical 
response, and the linearized hydrodynamical model then leads to coupled equations for the 
electrical field E and the current density J IT] [141 

V X V xE(r,t(j) = eothei (r,co)^E(r,(B) + /tOjUoJ(r,a)), (la) 
V [V • J(r, CO)] + J(r, to) = (j(r, co)E(r, co). (lb) 



(0 {(0 + i / Tomde) 



Also in linear response, the static electron density is homogeneous inside the metal while it 
is vanishing outside. In Eq. ([Tbl i. this electron density has been parameterized by the spa- 
tially piecewise constant AC conductivity a{r,(o) and TDmde is the Drude damping time that 
also occurs in the Drude local-response dielectric function £{(o) = 1 — (Op/[(o{(0 + i/ Tomde)], 
where ©p is the plasma frequency. We take parameters for gold, namely the plasma fre- 
quency h(Op = 8.812eV, Drude damping /j/Tomde — 0.0752eV, and Fermi velocity vp = 
1.39 X lO'^m/s, the same values as in Refs. ifTSl [T4l . In the Thomas-Fermi model and for 
ft'pTbmde ^ 1^ the nonlocal-response parameter (5 is proportional to the Fermi velocity through 
j3 = 1/3/(0 + 2) vp. Here D is the number of spatial dimensions from the point of view of the 
electron dynamics, which is the number of dimensions that are not quantum confined ll28l . For 
/3 — >^ we recover the local-response model where the dynamics is governed by Ohm's law with 
J = aE. We leave out the additional complexity of interband effects by taking eother(r, ft)) = 1 in 
Eq. ( [Tal l. The interband effects could be taken into account as well, following Refs. lfT2l[T3l[T5l . 

As stated in the Introduction, we do not consider the 'spill-out' of the electron density at the 
metal surface leading to quantum tunneling, as described by microscopic many-body calcula- 
tions 1251 129I . As an immediate consequence of this approximation, the normal component of 
the current J vanishes at the surface of the metal volume(s). We proceed to solve the coupled 
equations ([U self-consistently, with the usual Maxwell boundary conditons plus the additional 
boundary condition of the vanishing normal component of the current J at the metal surface(s). 
For further details of the model and the appropriate boundary conditions we refer to our recent 
theoretical work lfT4l . 

3. Computational method and implementation 

We first discuss the light extinction properties of infinitely long nanowire structures surrounded 
by free space. Since the electrons are confined in two directions but not quantum-confined in the 
nanowires, we should take D ~3 whereby the nonlocal parameter becomes ji — ^3/5vp 1281 . 
As the light source, we take a monochromatic plane wave incident normal to the wire axis and 
with an in-plane polarized electrical field (TM polarization). 

We solve Eq. ([T]i numerically with the aid of a finite-element method (FEM). To obtain a 
reliable and flexible implementation, we have built it as a nonlocal-response extension of a 
commercially available multi-purpose code for local-response (the RF module of Comsol Mul- 
tiphysics, version 4.1). The FEM is known to have a high ability to handle complex geometries, 
and to accurately model small surface details, such as gaps and tips ||30l . Furthermore, retarda- 
tion is automatically taken fully into account, which is important according to the very recent 
study by David and Garcia de Abajo fl5l . 

The translational invariance of the structure in one spatial dimension simplifies the calcula- 
tion, since the calculation domain (or 'grid') becomes two-dimensional. Within this plane, the 
wire geometry can be chosen at will. We imbed the 2D-projection of our metallic nanostruc- 
ture into a square computational domain. Perfectly-matched layers (PML) at the edges of this 




h(0 [eV] 

Fig. 1. Extinction-cross section CText (logarithmic scale) versus frequency for cylindrical 
nanorods for two radii: (a) radius a = 2nm, (b) a = 25 nm. (Text is normalized to the di- 
ameter of the rod. Both panels show comparisons of numerical simulations of Eq. l[T) to 
analytical results both for local response (/3 = 0) and for nonlocal response (jS = \/3/5vf). 
All numerical curves overlap the corresponding analytical curves. 

domain mimic the reflectionless coupling to the surrounding free space. For the meshing of the 
geometry we take advantage of the built-in algorithm of the software, paying special attention 
to mesh refinement needed to account for surface effects and abrupt changes in the surface 
topography. 

Our code runs on a pc. For details about our hardware implementation we refer to the Ap- 
pendix. The code was first successfully tested in the limit j3 ^ [see Eq. (fTbl)]. where it cor- 
rectly reproduces the local scattering response of various standard problems. Below, we report 
our benchmarking against analytical results, both for local and for nonlocal response. 

4. Benchmark problem: a single cylindrical nanowire 

First we compare our numerical method against analytical results for a single cylindrical 
nanowire, where the cylindrical symmetry allows for analytical solutions in terms of Bessel 
and Hankel functions, both for local and for nonlocal response lfT6l[T4l . Figure [T] summarizes 
the results of our benchmarking, and illustrates the strong dependence of the nonlocal effects 
on the subwavelength size of the nanostructure. The figure shows the dimensionless scaled 
extinction cross section (Text, defined as the cross section per length of the wire, divided by its 
diameter In more detail, Cext = (|^abs| + |^scat|)/(2a/o), where Iq represents the intensity of the 
incident plane wave, Pabs the absorbed and Pscat the scattered power per length of the wire, and 
2a is its diameter The powers are obtained by numerically integrating the Poynting vector on a 



circle surrounding the nanowire. 

Figure [Ha) compares the numerical results for a nanowire of radius a = 2nm to the exact 
analytical solution, Eq. (4) in Ref. 11411 . both for local and nonlocal response. For these tiny 
nanowires, nonlocal effect can be considerable. To give two examples, the relative difference 
of (Text for nonlocal against local response in the figure can be up to 15.33, and the resonant 
frequency is 6.20 eV for local and 6.39eV for nonlocal response, a considerable blueshift of 
A = 0.19eV, being more than 2 % of the resonance frequency. 

The analytical and numerical curves overlap almost completely for the local-response model, 
and likewise for the hydrodynamical nonlocal model. Hence only two of the four curves are vis- 
ible. Thus our numerical model accurately captures the prominent effects of nonlocal response, 
namely the blueshift of the surface-plasmon resonance «p/a/2 imillllTlini, as well as the 
confined bulk plasmon resonances above the plasma frequency Op lfT6l fT4l . More quantita- 
tively, the relative error of the numerically computed cross section for nonlocal response is 
always smaller than 0.4% in the entire frequency range of the figure, which includes many 
resonances, while for local response the relative error is always smaller than 0.6%. Further 
details about the accuracy and convergence of our numerical implementation are given in the 
Appendix. 

Figure[nb) shows the analogous four curves as Fig. [Ha), but now for a larger cylinder with 
radius a = 25 nm. As one can see, the four curves all overlap almost completely. Over the 
frequency range 1 .5 to 12 eV, the maximum relative difference in cross section for the nonlocal 
against the local model is 7%, which means that for most practical purposes nonlocal effects 
will be negligible for a single nanowire of this size. However, as we will see in Sec. |5] this 
conclusion does not carry over to a dimer of two such nanowires! The maximum relative error 
in the numerically computed cross section for nonlocal response in Fig. [TJb) is 0.89% while 
for local response it is 0.86%. These numbers show that the numerical method is accurate in 
a large frequency range also for larger nanostructures, here up to sizes where nonlocal effects 
can be neglected. 

Summarizing, the implementation for the cylindrical nanowire is accurate in a wide range 
of frequencies and length scales, correctly reproducing the location of resonances and their 
amplitudes even for high frequencies beyond the plasma frequency. 

5. Nonlocal effects in dimers of cylindical nanowires 

Having addressed single isolated nanowires, we now turn to nanowire dimers. For two closely 
separated nanostructures, the usual local-response model predicts a strong hybridization of the 
surface-plasmon resonance ll20l l22l . Thus dimers can even display resonances in the visible, 
with strongly enhanced fields in the tiny gap separating the two parts. Here we study the effects 
of nonlocal response both on the hybridization and on the field enhancement. We consider 
dimers of cylindrical nanowires (in this section) and bow-tie nanowires in Sec.|6l 

For definiteness, we consider a dimer of two identical cylindrical nanowires, with the same 
radius a = 25 nm as for the single nanowire of Fig. [Ttb), and separated by a few-nm gap of 
size d <^a. When excited by a plane wave normal to the wire axis and with the electrical field 
polarized across the dimer gap, the surface modes for this structure become highly confined 
to - and enhanced in - the narrow region of the gap. This field configuration is due to the 
hybridization of the modes of the isolated nanorods ll20ll21|[T2l . 

We examine the effects of nonlocal response on energy confinement and field enhancement 
for several gap sizes of 1 nm and beyond, where quantum tunneling, absent in our model, can 
be neglected 1251 . The results for cylindrical dimers are summarized in Fig.|2] We first discuss 
the three panels (al-a3) on the left of Fig.|2] that depict extinction cross sections for increasing 
gaps, before addressing the average field enhancements in the right-hand panels. 
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Fig. 2. Extinction-cross section versus frequency for a dimer of two cylindrical nanorods 
of radius a = 25 nm separated by a distance d, excited by a TM-polarized plane wave with 
wave vector perpendicular to the line connecting the centers of the cylinders, as illustrated 
in the inset with a/d not to scale. CText is normalized to the diameter of the single wire. The 
left panels labeled (a) depict (Text^ while the (b) panels on the right show the average field 
enhancement (y) as defined in the main text. The upper panels (al) and (bl) correspond to 
d = Inm, where the SPR appears at 3.28eV and 3.40eV for the local and nonlocal case, 
respectively. The middle panels (a2) and (b2) concern d = 2 nm where the SPR appears at 
3.63 eV in the local case and at 3.69 eV in the nonlocal one. The lower panels (a3) and (b3) 
correspond to d = 3nm. 



The results in panel (al) corresponding to a very small gap of c/ = 1 nm show a strong plas- 
monic interaction between the two nanowires. This interaction gives rise to a pronounced hy- 
bridization and consequently the main surface-plasmon resonance (SPR) now appears at much 
lower energy of 3.28 eV for local response and 3.40eV for nonlocal response, as opposed to the 
SPR in an isolated wire around 5.91 eV in both local and nonlocal case, recall Fig.lHb). The im- 
portant point is a pronounced nonlocal blueshift of 0.12eV of the hybridized dimer resonance. 
As the separation is doubled to d ~ 2nm, the hybridization decreases and consequently the 
SPR appears at higher frequency around 3.63 eV for local and 3.69 eV for nonlocal response. 
Again the effect of nonlocal response is a noticeable blueshift of the hybridized dimer reso- 
nance, still noticeable but smaller than for c/ = 1 nm. Finally, for d ^3 nm, the hybridization is 
again weaker, so that the SPR is again blueshifted in the direction of the single-nanowire SPR 
(see Fig. [I])- In all three cases, there is a nonlocal blueshift of the SP resonance frequency. This 
shift vanishes for larger d, because the blueshift for single nanowires with radius a = 25 nm 
also vanishes, as we saw in Fig.fljb). 

Now let us discuss the field enhancement in the right panels of Fig. |2] or more precisely the 
local field (intensity) enhancement factor /(r) = |E(r)p/|£'oP, whereE(r) is the local electrical 
field and £0 the amplitude of the incoming plane wave. Dimers support modes that are strongly 
localized in the gap separating the two nanowires. Due to the strong spatial localization, the 
amplitude of the local electrical field E(r) may by far exceed the amplitude Eq of the incoming 
plane wave. Rather than considering local field enhancements in single points, we will consider 
spatially averages, because local field probes such as atoms cannot be positioned with infinite 
precision. The average field enhancements (7) in Fig.|2b) were obtained by line-averaging /(r) 
over the narrow gap along the axis £ of the dimer, i.e. 



A direct comparison of left and right panels of Fig. |2] shows that the spectral dependence of 
the field enhancement (7) for lower frequencies is similar to the corresponding extinction cross 
section (Text- However, (7) peaks at higher frequencies than C7ext> both for local and nonlocal 
response. The agreement between the two types of curves is not complete, because the extinc- 
tion cross section is the sum of a scattering and an absorption cross section. The latter can be 
interpreted as a two-dimensional loss average inside the cylinders, whereas the field enhance- 
ment factor (7) of Eq. (I2]i is a more local one-dimensional spatial average of the empty space 
in between the cylinders. For that reason, the extinction aext can be high near 6eV while the 
field enhancement (7) is low. Indeed, there exists a resonant hybridized mode at this frequency 
(hence the peak in (Text) with a mode profile with low intensity on the line joining the cyUnder 
centers (which explains the low value for (7)); the mode intensity grows away from this line 
(mode profile not shown). 

For the narrowest gap of d=\ nm, Fig.|2|bl) shows a field enhancement of (7) = 1.4 x lO'* 
at the main SPR (i.e. at the lowest-energy peaks) for local response, and of (7) = 8.6 x 10^ for 
nonlocal response. Thus, the nonlocal response gives a strong average field enhancement, yet 
it is considerably smaller than for local response. As the gap distance increases to d = 2 nm, 
the maxima of (7) decrease for both types of material response. Finally in Fig|2b3) the gap 
has been increased further to d = 3 nm, and we still find considerable differences for nonlocal 
response, even though there is a further decrease both of the peak amplitudes of (7) and of their 
relative differences. 

Thus we find for the scattering off cylindrical nanowire dimers with radius a = 25 nm sepa- 
rated by 1 to 3 nm that nonlocal response effects are considerable, for cross sections but more 
so for field enhancements, even though in Sec.|4]we found that nonlocal effects were negligible 
for the single cylindrical nanowire with the same radius. 




(2) 



The above discussion compares nonlocal with local peak enhancements reduced amplitudes. 
Since the nonlocal peaks also shift in frequency, the nonlocal effects are even more pronounced 
when studying observables at fixed frequencies. For then Fig. |2|al) shows that (Text can be 
reduced by a factor of 3.1 at 4.0 eV, and enhanced by a factor of 3.4 at 3.7 eV. Nonlocal field 
enhancements for some frequencies can also turn out to be just a bit larger, for example by a 
factor 1.34 at 5.95 eV in Fig.[2jb3), even though we mostly find smaller nonlocal values for the 
field enhancement (7) in between the cylinders. Finally, it is interesting to notice frequencies 
in Fig. |2b) at which the local-response field enhancement peaks while the nonlocal-response 
field enhancement goes through a minimum. 

6. Nonlocal effects in bow-tie nanowires 

Let us now consider light scattering off the bow-tie nanowire dimer, the geometry as sketched 
in the inset of Fig. [3] Bow-tie structures can give rise to high field enhancement near the al- 
most touching sharp tips of the two triangles 131] [32l [33l |27l . Sharp features in nanoplasmonic 
structures are known to give strong field enhancements. The cylindrical dimer of the previous 
section did not have this type of field enhancement, and therefore it is interesting to compare 
nonlocal effects on field enhancement for the two types of dimers. 

In particular, we consider a dimer of two equilateral triangular nanowires with side L = 
45 nm. The tips of the triangles have been rounded with a radius of curvature of 1 nm. Such 
rounding for computational reasons is common practice ifTOl [34l [35]| . In Fig. [3] we present 
extinction cross sections and average field enhancements for bow-tie dimers, analogous to 
Fig. [2] The field enhancements are averaged over the line connecting the almost touching tips 
of the triangles. 

The results resemble those for the circular dimers, but with narrower and thus less overlap- 
ping plasmon resonances. One main feature is again that all (hybridized) surface-plasmon reso- 
nances are blueshifted due to nonlocal response. The case d ~l nm shows a strong interaction 
between the plasmons localized on the tip of the triangles, and the main SPR appears at 2.86 eV 
for local and 2.96 eV for nonlocal response. Larger distances imply weaker hybridization so that 
lowest-energy resonances shift to higher energies. Indeed for c/ = 2 nm the resonance frequency 
of the main SPR is 3.08 eV in the local model but blueshifted by 0.07 eV in the nonlocal one. 
Finally for d = 3nm the main SPR occurs at 3.24 eV in the local case and 3.29 eV in the non- 
local one. As for the dimer of cyhnders, for the bow-tie nanowire dimer we find larger nonlocal 
blueshifts of peaks in (Text in case of larger hybridization of surface-plasmon resonances. 

As to the importance of nonlocal effects in the field enhancement for the bow-tie dimer, 
Fig-EIb) shows that (7) peaks at 3.2 x 10^ in the local case and at 2.17 x 10"* in the nonlocal 
case for d ~ I nm. For the tip distance d ~ 2nm these values are 9.2 x 10^* and 7.2 x 10^, 
respectively, and they become 4.7 x 10^ and 4.6 x 10^ for d = 3nm. For field enhancement, like 
for cross sections, nonlocal effects turn out to be more important in case of larger hybridization. 

If instead of focusing on peak values we compare again nonlocal with local response at fixed 
frequencies, then for bow-tie dimers we can easily identify frequency intervals for which the 
field enhancement is larger for nonlocal than for local response, in contrast to what we found 
for the cylindrical dimers. For d = ?> nm at 4.81 eV, (7) is even 14.3 times larger for nonlocal 
than for local response. And for some frequencies, the nonlocal cross section is larger while for 
other frequencies it is smaller than the local response, and all differences are roughly within 
a factor of 5. Near 4.8 eV the field enhancements vanish even though the cross section peaks. 
Again this combination is a fingerprint of a resonant mode with low mode density near the 
center of the dimer 

The peak field enhancements (7) for the bow-tie dimers in Fig. [3b) are roughly a factor 2 
higher than for the cylindrical dimers in Fig.|2b), both for local and for nonlocal response. Thus 
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Fig. 3. Extinction-cross section versus frequency for a dimer of two equilateral triangles 
nanowires of side L = 45 nm separated by a distance d, excited by perpendicularly incident 
TM-polarized light. The inset shows a sketch with d/L not to scale. As in Fig. |2l the left 
and right panels show extinction cross sections CText and average field enhancements (y), 
respectively. Upper, middle, and lower panels again correspond to = 1,2 and 3nm. 

the quite different shapes of the dimers give rise to non-negligible but not too big differences in 
field enhancement, considering that all enhancement peaks (7) are of order 10"* or higher. This 
applies both to local and to nonlocal response. 



7. Discussion and conclusions 

We have implemented the hydrodynamical Drade model for arbitrary nanowire geometries as 
an extension of state-of-the-art numerical software in nanoplasmonics. Our code was tested 
against analytical results and was shown to be very accurate. We advocate the use of such 
benchmark problems, to be able to present results for more complex geometries with confi- 
dence. 

We studied tiny cylindrical nanowires, their dimers, as well as bow-tie dimers. In all cases we 
find that (hybridized) surface-plasmon resonances are blue-shifted due to nonlocal response. It 
is not simply the size of the plasmonic nanostructure that determines whether nonlocal effects 
are important. For example, we found that nonlocal effects were negligible for the extinction 
cross section of 25 nm cylinders, but important for their closely spaced dimers. 

We find that the usual explanations of plasmonic hybridization carry over to nonlocal the- 
ories. Below the plasma frequency we do not find new resonances due to nonlocal response, 
in agreement with our Ref. |[T4]|. but the hybridized resonances occur at higher frequencies 
than expected based on a local-response picture, with modified mode profiles. The nonlocal 
blueshifts are a correction to the larger hybridization energies. As an important conclusion, we 
find that the nonlocal blueshifts are larger for more strongly hybridized plasmonic structures. 
It is a package deal, so to say: one cannot have the one without the other, at least for dimers for 
which the individual parts are too big to exhibit any nonlocal blueshift. 

Bow-tie dimers have sharp tips and cylindrical dimers do not, but this did not result in large 
differences in nonlocal shifts in both cases. For the extinction cross sections of the dimers we 
found blueshifts in resonance peaks but hardly a change in their amplitudes. For the average 
field enhancements on the other hand, we find both blueshifts and a reduction in height of 
resonance peaks, roughly by a factor of 2 for cylindrical dimers and a factor of 1.5 for the 
bow-tie nanowires. Nevertheless this general conclusion is fully consistent with the fact that 
for some fixed frequencies, the field intensities in between the cylinders or bow-ties are extra 
enhanced due to nonlocal as compared to local response, by up to factors of 14. We also found 
frequencies for which the field enhancement peaks for nonlocal response but has minimum 
for local response. Thus it is important to take effects of nonlocal response into account in 
the context of spontaneous emission rates of nearby quantum dot emitters, fluorescence of dye 
molecules, and surface-enhanced Raman scattering (SERS) of bio-molecules, especially in the 
close vicinity of strongly hybridized plasmonic nanostructures. 

In this work we neglected quantum tunneling effects Il25ll29l . In their quantum many -body 
calculations, Zuloaga et al. |[25| identify a cross-over regime for dimer gaps between 0.5 and 
1 .Onm, where narrow -barrier quantum tunneling effects strongly reduce the classical hybridiza- 
tion energies, and a conductive regime for d < 0.5 nm. The classical (local-response) limit is 
also found in their calculations for large dimer separations. The message of our hydrodynami- 
cal calculations is that significant departures from classical local-response theory will already 
occur at larger dimer separations in the range 1 to lOnm, where quantum tunneling between the 
dimers is negligible, and that the local-response limit is found for large separations (i.e. for the 
individual a = 25 nm cylindrical wires for example). It would be gratifying to see experimen- 
tal nonlocal blueshifts for dimers, and to see many-body quantum calculations confirming not 
only the large-separation local-response limit, but also the hydrodynamical nonlocal-response 
blueshifts for dimers with separations in the range 1-3 nm as presented here. 
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9. Appendix: hardware implementation and accuracy studies 



All the presented results are obtained on a personal computer equipped with four Intel 2Ghz 
processors and 24Gb of RAM. Single-frequency calculations typically take about one minute 
for regular geometries with smooth boundaries. The frequency step used to make the spectra is 
0.0088 eV, and with this value all the data presented can be harvested in less than three days on 
our pc. 

The code that we have developed has undergone several reliability tests. In Fig.|4]we present 
a convergence study that has been carried out for the nanorod of radius a~2 nm with nonlocal 
response. The simulation box is a square of side W ~ 300 nm, surrounded by a PML that is 
30 nm thick. The walls of the box are in the far-field zone of the scatterer, so that they cannot 
influence the field scattered away from our system under test. 

First we studied the convergence of the nonlocally blueshifted value of the surface-plasmon 
resonance frequency, and we found that it converged to the analytical value 6.38 eV without 
much mesh refinement. Then at this resonance frequency, being a challenging spot, we studied 
the dimensionless convergence parameter 5 = |<7™t™ ~ '^ext\/'^e"t a function of the number 
of mesh elements used in the calculation. Meanwhile the number of mesh elements of both the 
box and the PML were kept fixed. The parameter 5 represents the relative error of <T™t™ with 
respect to the (7g"t, so that a vanishing 8 not just signifies convergence but rather convergence 
to the analytical value. 

Fig.Hdepicts 5 versus the number of mesh elements. It shows that the convergence sets in for 
less than 8000 mesh elements, a small number for modern pc's. In particular, 5 = 1.3% for 7241 
mesh elements and 0.12 % for 15106 mesh elements. This means that our code could easily run 
on a laptop, at least for the small structures considered here. The non-uniform mesh size allows 
also convergence on a pc for a larger cylinder with a — 25 nm, with 5 = 0.11% for 37830 mesh 
elements. Much larger structures can of course be handled on a powerful workstation. 
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Fig. 4. Relative error S of the numerically calculated extinction cross section cr"™' for 
a single cylinder of radius a = 2nm at the nonlocal-response resonant frequency tOsp = 
6.38 eV, and for a single cylinder of radius a = 25 nm at the frequency 5.91 eV, as a function 
of the number of mesh elements. Modeled with parameters of Au as given in Sec.|2] 



